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Purpose: In treatment planning of charged-particle radiotherapy, patient heterogeneity is conven- 
tionally modeled as variable-density water converted from CT images to best reproduce the stopping 
power, which may lead to inaccuracies in the handling of multiple scattering and nuclear interac- 
tions. Although similar conversions can be defined for these individual interactions, they would be 
valid only for specific CT systems and would require additional tasks for clinical application. This 
study aims to improve the practicality of the interaction-specific heterogeneity correction. 
Methods: We calculated the electron densities and effective densities for stopping power, multiple 
scattering, and nuclear interactions of protons and ions, using the standard elemental-composition 
data for body tissues to construct the invariant conversion functions. We also simulated a proton 
beam in a lung-like geometry and a carbon-ion beam in a prostate-like geometry to demonstrate 
the procedure and the effects of the interaction-specific heterogeneity correction. 
Results: Strong correlations were observed between the electron density and the respective effective 
densities, with which we formulated polyline conversion functions. Their effects amounted to 10% 
differences in multiple-scattering angle and nuclear-interaction mean free path for bones compared 
to those in the conventional heterogeneity correction. Although their realistic effect on patient dose 
distributions would be generally small, it could be at the level of a few percent when a carbon-ion 
beam traverses a large bone. 

Conclusions: The present conversion functions are invariant and may be incorporated in treat- 
ment planning systems with a common function relating CT number to electron density. This will 
enable improved beam dose calculation while minimizing initial-setup and quality-management of 
the user's specific system. 

Keywords: proton beam therapy, ion beam therapy, inhomogeneity corrections, multiple scattering, nuclear 
interactions 



I. INTRODUCTION 

In the present practice of radiotherapy treatment plan- 
ning, patient-specific material information is obtained 
from x-ray CT images. The CT number given to each 
image pixel represents the x-ray attenuation coefficient 
of the material of that pixel, from which the effective 
density for the treatment beam is estimated on the as- 
sumption of one-to-one correspondence. Taking water for 
a reference material, the effective density of a material is 
defined as the thickness ratio of water to the material 
for an equivalent dosimetric effect, which approximates 
to electron density for megavoltage photons or stopping- 
power ratio for charged particles. 

The empirical relationship between CT number and ef- 
fective density for a radiation of interest generally carries 
with it inaccuracies due to the limited tissue equivalency 
of sample materials. Schneider et al. developed the sto- 
ichiometric method to construct improved relations for 
body tissues^ They modeled the physics of x-ray attenu- 
ation using measured CT numbers of well-defined mate- 
rials, with which they calculated the x-ray attenuations 
in the standard body tissues of known elemental compo- 
sitions, mass density, and electron density^ in addition 
to their stopping-power ratios for protons, to construct 
the relations. 

With such conversion, patients are modeled as 



variable-density water to assure range accuracy for 
charged-particle beams in conventional heterogeneity- 
correction algorithms. As the body tissues deviate from 
water in elemental compositions, these algorithms only 
approximately address the other beam interactions. Mat- 
sufuji ct al. took the stoichiometric approach also for 
multiple scattering and nuclear interactions to study the 
validity of the conventional heterogeneity correction. 3 
They found that there would be 10%-level errors in 
multiple-scattering angle and nuclear-interaction mean 
free path of protons in bones although their effects on 
patient dose distributions would be small in realistic sit- 
uations. Nevertheless, these definite errors may be cor- 
rected by the interaction-specific density conversions. In 
fact, Szymanowski and Oelfkc proposed a specific CT- 
number conversion for multiple scattering^ Palmans and 
Verhaegen proposed another CT-number conversion for 
proton nuclear interactions. 5 

Those studies offer valid methods for improved hetero- 
geneity correction, which however have to be applied to 
individual scanning conditions of individual CT systems. 
The resultant conversion tables must be separately set 
up and selectively used in treatment planning systems 
with periodical constancy tests or updates, which would 
greatly complicate the system-management tasks in the 
clinical environment. The purpose of this study is to es- 
tablish a simpler procedure and to provide necessary data 
for equivalently accurate conversion. We take electron 
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density for a reference quantity and similarly construct 
conversions for the effective densities of body tissues as 
a function of electron density. As the conversion from 
CT number to electron density is commonly available in 
treatment planning systems, implementation of the in- 
variant conversion functions in algorithms will offer ac- 
curacy improvement without further complicating data 
management. 



main ingredient is water. As the u-dependent variation 
of the ps is within 1% under therapeutic conditions, 10 
we took the representative projectile speed v — 0.6 c or 
the nucleon kinetic energy E/A — 230 MeV to define the 
stopping effective density as projectile independent. 



C. Multiple scattering 



II. MATERIALS AND METHODS 
A. Electron density 

We used 92 of the ICRU body tissues,— excluding obso- 
lete, extreme, or artificially extracted materials such as 
ICRU-33 soft tissue, hydroxyapatite, calcifications, wa- 
ter, lipid, carbohydrate, cell nucleus, cholesterol, protein, 
and urinary stones. Electron density p c is calculated as 



(1) 



where Zi and A Ti are the atomic number and the atomic 
weight of element i, (Z/A r ) w = 0.5551 is the mean Z/A r 
of water, and Wi and /9 m are the elemental mass fraction 
and the mass density of the material. We represent p m , 
Pc and all other densities as non-dimensional ratios to 
those of water in this study. 



B. Stopping power 

The Bethe theory leads to the stopping-power ratio of 
a material to water j£ or the stopping effective density, as 
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where m = 0.511 MeV/c 2 is the electron mass, v and 
c are the speeds of the projectile and light, and I and 
I w are the mean excitation energies of the material and 
water. We calculated the I values of the body tissues 
with the Bragg rule 



■A Ti 



(3) 



where elemental Jj/eV values were H:19.2, C:81, N:82, 
0:106, F:112, Na:168, Mg:176, P:195, Cl:180, K:215, 
Ca:216, and Fe:323 for constituents of solid or liquid 
compounds^ This leads to I w — 75.3 eV for water. Al- 
though recent studies suggest slightly higher I w values,— 
we used this J w value in the denominator of Eq. ^ to 
cancel out the common systematic error of the compound 
/ value in the numerator for general body tissues whose 



Scattering power T is the increase of Gaussian angular 
variance per transport distance.- In some formulations, 
it is in the form of 
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where E s — 15.0 MeV is a constant energy, p is the par- 
ticle momentum, and X is the scattering length of the 
material. Gottschalk derived the scattering length for 
heavy charged particles?^ For a composite material, it is 
given by 
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This leads to X w — 46.88 cm for water. We define the 
scattering effective density as the scattering-power ratio 
of the material to water, 



PT T w X 
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D. Nuclear interactions 



Sihver et al. made an empirical modification to 
the geometric model for nucleus-nucleus collision cross 
section; 12 ! 13 resulting in 
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for ions, 



(8) 



where we consistently use symbols A for mass number 
of the projectile and A r for atomic weight of the target 
and introduced tq = 1.36 fm for effective nucleon radius 
and bo for overlap parameter. In their model, the energy 
dependence of the cross section would be almost common 
among target body-tissue materials in the therapeutic 
energy region and is thus ignored in this study. 

For a compound material of mass density p m , the num- 
ber of elemcnt-i nuclei per volume is proportional to 
PmWi/A ri . The total nuclear cross section per volume 
is thus proportional to J2i a NiPmWi/ A ril which is pro- 
portional to the incidence of nuclear interactions. We 
define the nuclear effective density as the ratio of the in- 
cidence of nuclear interactions in the material to that in 
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water, resulting in 
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where onh an d ^no are the cross sections of hydrogen 



and oxygen nuclei and M T - 



18.015 is the molecular 



weight of water. Unlike the other effective densities, we 
formulated the nuclear effective density as projectile de- 
pendent. 



TABLE I: Polyline conversion functions from electron density 
p c to stopping effective density ps, scattering effective density 
pT, and nuclear effective density pN for protons and ions of 
4 < A < 16. 



Pc 







0.9 0.9 1.035 1.4 2.0 



ps/p a 1.000 1.000 1.028 1.000 0.973 0.942 

Pt/pc 0.995 0.995 0.77 0.995 1.32 1.66 

Pn/Pc for protons 0.992 0.992 1.07 0.992 0.89 0.78 

Pn/Pc for ions 0.987 0.987 1.12 0.987 0.81 0.64 



III. RESULTS 

E. Example of corrective conversions 



To demonstrate the procedure and to evaluate the 
significance of the interaction-specific corrective conver- 
sion, we carried out simplistic proton and carbon-ion 
beam-transport calculations in one-dimensional geome- 
tries, each with and without the corrective conversions. 
We applied a proton beam of range 9.8 cm in water to 
a lung-like geometry consisting of staked slabs of 1-cm 
soft tissue, 2-cm bone, 1-cm soft tissue, 8-cm lung, and 
3-cm tumor. We applied a carbon-ion beam of range 
22.4 cm in water to a prostate- like geometry consisting 
of stacked slabs of 6-cm soft tissue, 6-cm bone, 4-cm soft 
tissue, and 4-cm tumor. We assigned electron densities 
1.0 for soft tissue and tumor, 1.4 for bone, and 0.25 for 
lung and set the beam ranges to coincide with the tumor 
depths without the correction. These beams were point- 
like mono-directional at the phantom entrance and were 
transported through the slabs at 1-mm intervals with 
quantities of interest, 



d(x) = / p(x')dx' (11) 
Jo 




evaluated as a function of transport distance x, where 
generic depth, beam size, and density, (d,a y ,p), repre- 
sent the effective depths and the beam sizes calculated 
with the respective effective densities. For the integra- 
tion, the Fermi- Eyges theory in the numerical formic was 
implemented with a semi-empirical scattering power for 
water, T W J^ We are interested in the beam size corrected 
for multiple scattering, <7 yT , and the effective depth cor- 
rected for nuclear interactions, g?n, to compare to the 
effective depth and beam size in the conventional het- 
erogeneity correction, (ds,<J ys )- The effective depths are 
geometric parameters that will not depend on beam prop- 
erties except for the c?n that is slightly different between 
protons and ions. We specifically refer to the effective 
depths of the distal end of the target as the effective 
ranges. 



A. Corrective conversion functions 



Figure [T] shows correspondences between the elec- 
tron density and the effective densities of body tis- 
sues for proton and ion beams. There was a con- 
centration of tissues around (p c , ps/pe, PT/Pe, Pn/Pc) = 
(1.035,1.000,0.995,0.99). The low density (p e = 0.258) 
for the lung tissue was attributed to the air content. 
Generally in the p c > 1 region, a weak negative corre- 
lation between ps/p e and p e was observed. That was 
attributed to the low / values of carbon-rich adipose 
tissues in the low p c region and the high / values of 
calcium-rich bone tissues in the high p e region. Sim- 
plistically for a single- element material, Coulomb scat- 
tering depends on p m Z 2 /A T whereas the electron density 
depends on p m Z/A r . As there is approximate linear rela- 
tionship among p m , A T , Z, and p c for common materials, 
the ratio pt/ p e oc Z showed a strong positive correlation 
with p c . Similarly, because the nuclear effective density 
depends approximately on p m A T ~ 1 ' 3 for large A r , the ra- 
tio pn/Pc oc A r 2 / 3 /Z showed a negative correlation with 
p c . The pn/pe ratios for protons largely deviated from 
those for the ions due to the distinctive formulation of 
the overlap parameter bo in Eq. ([5]). 

For the conversion functions, we set a discontinuity 
point at p c — 0.9, where none of these tissues existed, 
an inflection point at the center of the tissue concen- 
tration, and another inflection point at p c = 1.4 in the 
bone-tissue region. Similarly to the lung tissues, tissue- 
air mixing may occur on the surface of skin and epithelial 
tissues due to finite spatial resolution of the CT image, 
for which we assigned constant p/p e ratios in the p c < 0.9 
region. Although material mixing may generally confuse 
any conversion schemes, its influence is limited as there 
normally are only a few tissue boundaries along the beam 
path. Because the variation of pn/Pc among the ions was 
small, we defined one conversion function for the ions of 
4 < A < 16. Table U shows the resultant conversion fac- 
tors as polyline functions, which are also shown in Fig.[T] 
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FIG. 1: Ratios of (a) stopping effective density ps, (b) scat- 
tering effective density pT, (c) nuclear effective density pN for 
protons, and (d) nuclear effective density pn for helium ions 
(x), carbon ions (A), and oxygen ions (+) to electron density 
p e with the polyline conversion functions (solid lines). 



FIG. 2: Corrections of stopping and nuclear effective depths, 
(ds — d c , g(n — d c ), and uncorrected, conventional, and cor- 
rected beam sizes, {ct Vc , o vs , c yT ), with background colors 
indicating electron densities, (a, b) for a proton beam in a 
lung-like geometry and (c, d) for a carbon-ion beam in a 
prostate-like geometry. 



B. Example of corrective conversions 

Figures EJa) and [2fc) show the development of the 
depth corrections for stopping and nuclear interactions, 
(ds — d c ) and (c?n ~ d e ), for the proton beam in the 
lung-like geometry and for the carbon-ion beam in the 
prostate-like geometry. The corrections of the stopping 
effective depth, which would assure the range accuracy, 
were within about ±0.1 cm in both cases. The bone 
thickness mainly determined the corrections of the effec- 
tive ranges. The nuclear effective range (Rn = 21.6 cm) 
deviated from the stopping effective range (i?s = 22.3 
cm) by —0.7 cm in the prostate-like geometry. 

Similarly, figures HJb) and [UJd) show the development 
of the beam sizes. The uncorrected beam size a y and the 
conventional beam size a vs were almost indistinguishable 
in these cases, while the corrected beam size a VT devi- 
ated from them by about 5% for the proton beam in the 
lung-like geometry. The bone thickness to scatter and 
the distance to travel determined the amount of the cor- 
rection. 



IV. DISCUSSION 

A distinctive feature of this study is the generality of 
the conversion functions. For example, Szymanowski and 
Oelfke constructed a conversion function for the angular 
scaling factor or \J Pt/ps m our notation. 4 While their 



function was only valid for a specific CT system, the p c - 
Pt function in Table|T]is generally valid. As the biological 
data and the physical models are common or equivalent, 
the direct CT-number conversion and the corrective con- 
version through p G are essentially equivalent in accuracy, 
assuming that the corrective conversion is used with a 
stoichiomctrically constructed CT-number-p e tablei^ 
Intrinsic error of CT number, such as that from x-ray 
hardening effect (typically 1% for soft tissues and 2% for 
bones), 16 influences both conversion schemes in a similar 
way. 

Matsufuji et al. found that, for a typical bone of p c « 
1.35, the errors of the conventional heterogeneity cor- 
rection would be —10% in multiple-scattering angle and 
— 10% in mean free path of proton-nuclear interaction. 3 
For a tissue of p c — 1.35, the polyline functions in ta- 
ble U give pt/Pc = 1.21 and p^/pc — 0.904 for pro- 
tons. Because the multiple-scattering angle is propor- 
tional to y/pr and the mean free path is proportional to 
1/pN, the pt and pn conversions would make +10% and 
+11% corrections to multiple-scattering angle and mean 
free path, respectively. The +10% and +11% corrections 
correspond to 9 % and 10% effects in the corrected val- 
ues, which agree with their original estimations. They 
also simulated realistic situations, in which bone was not 
dominant in the beam path, and concluded that the ef- 
fects on patient dose distributions would be small, which 
we confirmed in our examples. 

The nuclear interactions cause attenuation of primary 
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particles. Lee et al. proposed an approximation for- 
mula for proton fluence <$>(R) = $(0) (1 + 0.012i?/cm) 
as a linear function of residual range R in water,— to 
which the nuclear effective range i?N should be applied 
in the interaction-specific heterogeneity correction. For 
the prostate-like geometry ignoring the small i?N differ- 
ence between protons and ions, the error of the nuclear 
effective range with the conventional heterogeneity cor- 
rection (AR = i?s — Rn = 0.7 cm) would cause a rel- 
ative error of —0.7% to the proton survival $(0)/$(i?) 
and consequently to the dose to the tumor. Palmans 
and Verhaegen reported up to 0.5% dosimetric effects for 
their examples that contained less bone than ours£ For 
carbon-ion beams, the dosimetric effect would be a few 
or more times larger for their large nuclear cross section, 
i.e. a few percent, which could be clinically influential. 

The fragmentation processes in ion beams are further 
complex. In the participant-spectator model, the target 
nucleus will not directly influence how the projectile nu- 
cleus may break up J£ There are some experimental data 
for the relative yields of projectile fragments, 19 which 
may be thus reasonably invariant to target materials. 
Due to experimental difficulties in precise measurement 
of a treatment-beam spectrum, Monte Carlo simulation 
will be useful for building a detailed beam model in nu- 
merical or analytical form.— 



V. CONCLUSIONS 

For general body tissues, we found strong correlations 
between electron density and the effective densities that 



characterize the strengths of stopping power, multiple 
scattering, and nuclear interactions of the projectile pro- 
tons and ions. The stopping effective density deviated 
from the electron density by up to a few percent and the 
scattering and nuclear effective densities deviated by up 
to a few tens percent, which are consistent with other 
studies. To correct those errors, invariant corrective con- 
version functions from the electron density to the indi- 
vidual effective densities were constructed. Their effect 
on patient dose may be negligibly small in general, al- 
though it could be at the level of a few percent when an 
ion beam traverses a large bone. 



The electron density should be derived from the CT 
number with a conversion table that has been con- 
structed with the best accuracy, i.e. stoichiometrically, 
for treatment planning of proton and ion-beam therapy. 
The proposed extension of heterogeneity correction will 
enable improved beam dose calculation without seriously 
sacrificing simplicity or efficiency of the patient model 
and the algorithms, while quality-management tasks for 
CT and treatment panning systems will be kept as they 
are in the current clinical practice. 
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